Silica-Like Sequence of Anomalies in Core-Softened Systems 
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In this paper we present a simulation study of density, structural and diffusion anomalies in core- 
softened system introduced in our previous publications. It is well-known, that with appropriate 
parametrization, core-softened systems are remarkable model liquids that exhibit anomalous proper- 
ties observed in tetrahedral liquids such as silica and water. It is widely believed that core-softened 
potentials demonstrate the water-like sequence of anomalies. We show that with increasing the 
depth of the attractive part of the potential the order of the region of anomalous diffusion and the 
regions of density and structural anomalies is inverted and have the silica-like sequence. We also 
show that the slope of the Widom line is negative like in water. 

PACS numbers: 61.20.Gy, 61.20.Ne, 64.60.Kw 
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In the last decades there is a growing interest to a 
class of fluids that display anomalous thermodynamic 
and kinetic behaviors. The most common and well known 
example is water. The water phase diagram has re- 
gions where a thermal expansion coefficient is negative 
(density anomaly), self-diffusivity increases upon com- 
pression (diffusion anomaly), and the structural order 
of the system decreases with increasing pressure (struc- 
tural anomaly) [3-@|. Later on it was discovered that 
many other substances also demonstrate similar behav- 
ior. Some typical examples are silica, silicon, phosphorus 
and many others 

As it was found in experiments [l| and simulations 
0, 0] , the water anomalies have a well-defined sequence: 
the regions where these anomalies take place form nested 
domains in the density-temperature |3| (or pressure- 
temperature Q) planes: the density anomaly region 
is inside the diffusion anomaly domain, and both of 
these anomalous regions are inside a broader structurally 
anomalous region. This water-like behavior was found in 
systems with spherically symmetric core-softening poten- 
tials with two length scales ll|427|. 



However, in other anomalous systems the sequence of 
anomalies may be different. For example, in computer 
simulation of the system with the Van Beest-Kramer-Van 
Santen (BKS) potential the hierarchy of anomalies for 
silica is different compared to water [7[. In this case the 
diffusion anomaly region contains the structural anoma- 
lous region which, in turn, has the density anomaly re- 
gion inside. To our knowledge this is the only example 
of such inversion of the order of the anomalies discussed 
in literature till now. 

It is widely believed that in the core-softened systems 
hierarchy of anomalies has the water-like type. For exam- 
ple, Yan et al. [13, EH characterized the structural, ther- 
modynamic, and kinetic properties of a family of discon- 



tinuous core-softened potentials that vary in the length 
scale of the soft repulsion region. In each case, they found 
the same relationship between anomalous regions as ob- 
served for water: structural anomalies preceded diffusiv- 
ity anomalies, which preceded density anomalies. In Ref. 
[19j it was discussed an approach for the analysis of the 
anomalous behavior based on the well known Rosenfeld 
scaling relations |28[ which connect the transport coef- 
ficients with excess entropy. In this case there was ob- 
tained the explicit equation determining the appearance 
of the anomalies in system [19(: 
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where p is the number density, Sex = S — Sid is the 
excess entropy, equal to difference between total S and 
ideal gas Sid entropies, and c is the property-specific con- 
stant. Based on Rosenfeld's scaling parameters [19], it 
was shown that value c = corresponds to structural 
anomaly, c = 0.42 to diffusion anomaly, and c = 1 - 
to density anomaly. From ([T]) one can conclude that 
anomalous behavior always occur in the water-like order: 
structural anomalies precede diffusivity anomalies, which 
in turn precede density anomalies. This conclusion was 
verified by the computer simulation of the system inter- 
acting through a two-scale potential introduced by Jagla 
[29 I. However, as it was shown in our previous publi- 
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Rosenfeld scaling fails in the vicinity of 
anomalies, and it can not be used for the correct analysis 
of the order of the anomalies. 

It is important to note that the equation ([TJ contains 
two conditions which are basing on strict thermodynamic 
arguments: the conditions for structural and density 
anomalies (l9| . This means that the density anomaly 
should always follow after the structural anomaly. On 
the other hand, the diffusion anomaly can be located at 
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TABLE I: The potential parameters used in simulations (Eq. 
©)• 



every place. For example, for silica the diffusion anomaly 
precedes the structural and density anomalies. 

This paper presents a simulation study of anomalies 
in core-softened system introduced in our previous pub- 
lications (20I - I27I ] . We investigate the sequence of the 
anomalous regions and find that with increasing the at- 
tractive part of the potential the system demonstrates 
both water-like and silica-like behavior. It is also shown 
that the potential qualitatively correctly reproduces the 
behavior of the Widom line for the liquid-liquid transi- 
tion. 

In the present study we investigate a system of parti- 
cles interacting via the potential with "hard" core, repul- 
sive shoulder and attractive well (21. 25|. 

The general form of the potential is written as 



U(r) 
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Here k\ — k2 — 10.0 and the parameters of the poten- 
tials ([2]) are given in Table 1 . The family of the potentials 
with <7i = 1.35 and different attractive wells is shown in 
Fig. Q] It should be noted that the potential ([2]) is very 
similar to the Fermi- Jagla potential suggested recently in 
Ref. (32j . where the Fermi distribution function is used 
except for the hyperbolic tangent in Eq. ([2|) in order to 
describe the smoothed step [33 1. 

In the remainder of this paper we use the dimensionless 
quantities: f = r/a, P = Pa 3 /e, V = V/Na 3 = l/p,f = 
ksT/e. As we will only use these reduced variables, we 
omit the tildes. 

In Refs. it was shown that these systems 

demonstrate anomalous behavior. A relation between 
phase diagram and anomalous regions was also discussed 
in these articles. Our later publications gave detailed 
study of diffusion, density and structural anomalies in 
this system (2H - E7I] . 

In the present article we carry out a molecular dynam- 
ics and Monte Carlo study of the core-softened systems 
and monitor the change in phase diagram and anomalous 
regions with changing the potential parameters. The de- 
tails of simulations can be found in Ref. |2f| . 

In order to find the melting lines we carry out the free 
energy calculations for different phases and construct a 
common tangent to them. For our potentials we com- 
puted the free energy of the liquid by integrating the 
equation of state along an isotherm J34j]: F ^' fc F ^ ( - p - ) = 
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FIG. 1: (Color online) Family of the potentials with o\ = 1.35 
and different attractive wells. The curves are numerated in 
accordance with Table 1. 
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tal phases were determined by the Monte Carlo simula- 
tions with the method of coupling to the Einstein crystal 
34]. In this case the excess entropy can be computed 
via S ex — ^vfc F j 1 where U is the internal energy (34j . 
The total entropy is S = S ex + Sid, where the ideal gas 
|ln(T)~ln(p)+ln(^™^) 3/2 - 5/2 ' 



entropy is 



h 3 



'-)■ 



Here we present the anomalous regions for the three 
systems shown in Fig. [TJ 

The density anomaly means that density increases 
upon heating or that the thermal expansion coefficient 
becomes negative. Using the thermodynamic relation 
{DP / dT) v = ckp/Kt, where ap is a thermal expansion 
coefficient and Kt is the isothermal compressibility and 
taking into account that Kt is always positive and fi- 
nite for systems in equilibrium not at a critical point, we 
conclude that density anomaly corresponds to minimum 
of the pressure dependence on temperature along an iso- 
chor. This is the most convenient indicator of density 
anomaly in computer simulation. 

Initially the structural anomaly was introduced via or- 
derparameters characterizing the local order in liquid 
S H El Ei-El]. However, later on the local order 
was also related to excess entropy of the liquid which 
is defined as the difference between the entropy and the 
ideal gas entropy at the same (p, T) point: S ex = S — Sid- 
In normal liquid excess entropy is monotonically decaying 
function of density along an isotherm while in anomalous 
liquids it demonstrates increasing in some region. This 
allows to define the boundaries of structural anomaly at 
given temperature as minimum and maximum of excess 
entropy. 

The behavior of the diffusion coefficient, pressure and 
excess entropy for the purely repulsive system with o± = 
1.35 has been discussed, for example, in Refs. 2(| 2^, 25h - 
271 ] . One can see that all three anomalies take place in 
the system. Fig. [U[a) places the regions of the anomalies 
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in the phase diagram. As it was shown in [20, 21. 26| 
(see, for example, Fig. 1 in Ref. [2l]|). the phase diagram 
of the system consists of the high density and low den- 
sity Face Centered Cubic (FCC) phases, corresponding 
to hard core and repulsive shoulder parts of the poten- 
tial, separated by a sequence of crystalline phases. In 
Fig. [2ja) we show the low density part of the phase dia- 
gram with the FCC and FCT (Face Centered Tetragonal) 
phases. 

One can see that the anomalous regions correspond 
to the picture proposed for water i.e. the diffusion 
anomaly region is inside the structural anomaly and the 
density anomaly is mainly inside the diffusion anomaly. 

In Ref. 26[ it was shown that in the system with 



the purely repulsive potential the diffusion and density 
anomalies inverted with respect to each other with in- 
creasing the repulsive core diameter, i.e. now diffusion 
anomaly region is inside the density anomaly one. 

Next we consider the influence of attraction on the 
anomalous behavior of the system. For this, we study the 
system with step size a\ = 1.35 and different well depths 
(Table 1 and Fig. [1]). In this case one can see that with 
increasing the depth of the attractive well the sequence of 
the anomalies also inverted: for w = we have the water- 
like order of the anomalies (see Fig. [2a)), at w = 0.2 the 
locations of the diffusion and structural anomalies almost 
coincide (Fig. d[b)). At last, for w = 0.4 the diffusion 
anomaly region contains the structural anomalous region 
which, in turn, has the density anomaly region inside 
(Fig. Hfc)) 40j. As a result, for w = 0.4 we obtain the 
configuration of anomalies which is the same as in silica 

It seems that a most popular point of view is that the 
hypothesized liquid-liquid critical point is the thermo- 
dynamic source of all water anomalies 3l|, I4ll443l | . the 



terminal point of a line of first-order liquid-liquid phase 
transition. The line emanating from this critical point 
is sometimes called the Widom line and is often consid- 
ered as an extension of the coexistence line into the one- 
phase region (see, for example, 3l|, [32J, |45|).This line is 



determined by the lines of the maxima of the thermody- 
namic response functions which asymptotically approach 
one another as the critical point is approached 3l[ |45[ . 
The lines of the liquid-liquid rfiase transition in the com- 
puter simulations of water 3l| and silicon 44j have the 
negative slope. On the other hand, earlier computer sim- 
ulations of the isotropic core-softened potentials (see, for 
example, [3l|) suggested that the slope of the Widom line 
is positive in these systems. 

In Fig. ([3]) we show the maxima lines of isobaric heat 
capacity Cp along isobars and isotherms, the line of the 
isothermal compressibility maxima Kt, the line of the 
thermal expansion coefficient ap minima along with the 
melting line for the systems 2 and 3 (Table 1). One can 
see, that for our potential the Widom line has a slope 
which coincides with the results for water |31j and silicon 
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FIG. 2: (Color online) (a) Location of anomalous regions at 
the low density part of the phase diagram of the system with 
<7i = 1.35, where Face Centered Cubic (FCC) and Face Cen- 
tered Tetragonal (FCT) phases are shown. The water-like 
order of anomalies takes place: the density anomaly region 
is inside the diffusion anomaly domain, and both of these 
anomalous regions are inside a broader structurally anoma- 
lous region. Location of anomalous regions at the phase dia- 
gram for (b) system with o\ = 1.35 and w — 0.2; (c) ai = 1.35 
and w = 0.4 with the silica-like order of anomalies (see Table 
I)- 
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FIG. 3: (Color online) The maxima lines of isobaric heat ca- 
pacity Cp along isobars and isotherms, the line of the isother- 
mal compressibility maxima Kt, the line of the thermal ex- 
pansion coefficient minima along with the melting line for the 
systems 2 (a) and 3 (b) (Table 1). 



44( 1 . It should be noted that the similar results were 
obtained recently [32[ for the potential similar to ([2]). 

In conclusion, this publication represents a detailed 
computer simulation study of anomalous behavior of 
core-softened systems proposed in our publications (20I — 
271 ] . The core-softened potentials are widely investigated 
because they reproduce the water-like anomalies. Taking 
into account that the anomalies also exist in the systems 
where the hydrogen bonds are absent, it seems that the 
unusual properties of water are quite universal and can be 
investigated with the help of the isotropic core-softened 
potentials. To our knowledge, all isotropic core-softened 
potentials which were considered in the previous publi- 
cations show the sequence of anomalies characteristic for 
water: the density anomaly region is inside the diffusion 
anomaly domain, and both of these anomalous regions 
are inside a broader structurally anomalous region. On 
the other hand, in other anomalous systems the sequence 
of anomalies may be different. For example, the hierarchy 
of anomalies for silica is different compared to water. In 
the present article we analyze the possibility of changing 
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FIG. 4: (Color online) Schematic of the regions wherein struc- 
tural, diffusivity, and density anomalies are found within the 
temperature-density planes. Region of anomalous density be- 
havior always appears as nested dome within the structural 
anomalous envelope while the diffusivity anomaly may be (a) 
between the structural and density ones (see Fig. [5] (water- 
like behavior); (b) inside the density anomalous region; and 
(c) the outermost envelope (silica- like behavior). 



the order of anomalies regions depending on the param- 
eters of the potential. It is shown that for the potential 
@ for small values of the repulsive step the sequence 
of anomalies is the same as in water, however, with in- 
creasing the width of the repulsive shoulder the order of 
the region of anomalous diffusion and the region of den- 
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sity anomaly is inverted. With increasing the depth of 
the attractive well we obtain the configuration of anoma- 
lies which is the same as in silica Q , where the diffusion 
anomaly region contains the structural anomalous region 
which, in turn, has the density anomaly region inside. It 
seems that this is the first case when the isotropic core- 
softened system demonstrates the sequence of anomalies 
which is different from the water one and may be the 
same as in silica. It is also shown that for the potential 
<j2j) the slope of the Widom line is negative as in the case 
of water and silicon. 
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